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Abstract 



We extend the effective theory of Minimal Lepton Flavor Violation (MLFV) 
by including four-lepton operators. We compute the rates for fi — > 3e and r — > 3£ 
decays and point out several new ways to test the hypothesis of MLFV. We also 
investigate to what extent it will be possible from (future) experimental informa- 
tion to pin down the contributions of different effective operators. In particular we 
Qh! look for experimental handles on quark-lepton operators of the type Ifilj x qTq by 

Q_ ( ' working out their contribution to hadronic processes such as r —* /U7r°, 7r° — * fie, 

T — > tu, as well as to purely leptonic decays such as u — > 3e through loop effects. 



1 Introduction 



In a recent work |Tj we have extended the notion of Minimal Flavor Violation (MFV) j2HHl 
H] to the lepton sector of beyond the Standard Model (SM) theories. The MFV hypothesis 
states that the irreducible sources of (lepton) flavor symmetry breaking are linked in 
a minimal way to the structures generating the observed pattern of fermion masses 
and mixing. While this idea has a straightforward and unique realization in the quark 
sector 01313] (the SM Yukawas are the only sources of quark-flavor symmetry breaking), 
the situation in the lepton sector is different, mainly due to the possibility /necessity to 
break the U(l) symmetry associated with total lepton number. 

The MFV framework has two particularly attractive features. On one hand, it implies 
a suppression of FCNC processes induced by new degrees of freedom at the TeV scale to 
a level which is consistent with present experimental constraints. Moreover, it provides a 
predictive framework that links the FCNC couplings to the fermion spectrum and mixing 
structure (up to an overall scale factor). For leptons it relates lepton- flavor mixing in 
the neutrino sector to lepton-flavor violation in the charged lepton sector. 
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In ref. |T] we found two ways to define the sources of flavor symmetry breaking in the 
lepton sector in a minimal and thus very predictive way: i) a scenario where the left- 
handed Majorana mass matrix is the only irreducible source of flavor symmetry breaking; 
and ii) a scenario with heavy right-handed neutrinos, where the Yukawa couplings define 
the irreducible sources of flavor symmetry breaking and the right-handed Majorana mass 
matrix has a trivial flavor structure. 

For each scenario, according to the MLFV symmetry principle, we have constructed 
the basis of dimension six operators contributing to processes with only two charged 
leptons |L. In this context we have investigated the sensitivity of processes such as 
£% £jl and // — > e conversion to the scales of lepton flavor (Alfv) and lepton number 
(Aln) violation and we have pointed out distinctive predictions of MLFV as a function 
of S13 and the CP violating phase S. 

In this paper we extend previous work in several respects: 

• we enlarge the MLFV effective theory to include operators contributing to four- 
lepton processes and we give predictions for decays such as /i — > 3e and r — > 3£. 

• we investigate to what extent it will be possible to reconstruct from (future) ex- 
perimental information the dynamics of a given MLFV model, which amounts to 
pin down the relative size of the couplings appearing in the effective lagrangian. 

We tackle this second issue by working out a number of predictions for LFV decay 
rates. In particular, we study how LFV decays involving hadrons can be used to probe 
operators involving two leptons and two quarks (such as £iT£j x qTq). We then show 
how, under suitable assumptions, the decay fi — > 3e can also be used to probe operators 
of the type £iT£j x qTq through loop effects. 

The paper is organized as follows. In section|2]we recall the two realizations of MLFV 
introduced in Ref. pQ and give the complete basis of dimension six effective operators, 
including four-lepton operators. In section |3] we present results for the rates of various 
£ — > £'£"£"' decays at tree level in the effective field theory, explore their phenomenology, 
and point out new testable predictions of the MLFV scenario. In section 0] we discuss 
LFV decays involving hadrons, while in section El we calculate the loop-induced contri- 
bution to fi — > 3e from four-fermion operators involving two quark fields. We present 
our conclusions in section |H1 Details concerning the phase space integration, the renor- 
malization group analysis of section |5] and the matching to chiral perturbation theory 
are given respectively in appendices IAIBI and 

2 Extended Operator Basis 

In this section we recall the two scenarios of MLFV identified in Ref. Q and present 
the dimension six operator basis, extended to cover four-lepton processes. In order to 
formulate the minimal flavor violation hypothesis for leptons, one needs to identify the 
field content of the theory, the flavor symmetry group and the irreducible sources of 
symmetry breaking. Two possibilities arise (see ref. £Q for details). 

Minimal field content: In this scenario the flavor symmetry group is Glf = SU(3)l x 
SU(3)e, acting on three left-handed lepton doublets L l L and three right-handed charged 
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lepton singlets e R (SM field content). The breaking of the C/(1)ln is independent from 
the breaking of Glf and is associated to a very high scale A L n- The irreducible sources 
of lepton-flavor symmetry breaking are X l J and g l J , defined by 1 

4jym.Br. = ~K j eUH^Li) - -J- g^{Lir 2 H){H T T 2 Li) + h.c. (1) 

^LN 

Extended field content: In this scenario the maximal flavor group is Gxf x SU(3) Ur , 
acting on three right-handed neutrinos, u l R , in addition to the SM fields. The right- 
handed neutrino mass term is flavor diagonal (M' l J = M u 5 tj , \M V \ 3> v) and its effect 
is to break £/(1)ln as well as to break SU(3) UR to 0(3) UR . The remaining lepton-flavor 
symmetry is broken only by two irreducible sources, \ l J and \ % J, defined by 

4jym.Br. = -A« e R (H^L{) + i^V R (H T r 2 L{) + h.c. (2) 

In most SM extensions, at some scale Alfv above the electroweak scale and well 
below Aln (or M v ) there are new degrees of freedom carrying lepton flavor quantum 
numbers. As long as the underlying model respects MLFV, at scales below A L fv the 
physics of lepton flavor violation is described by the following effective lagrangian (up to 
higher dimensional operators) 
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A 2 

LFV i=1 * ' LFV \ j=1 



with operators defined by: 



0« = g'H^e R a^\ e AL L B„ v 



= l l yal l HUD^H 



Of L = L Ll »AL L QlIuQl 



0^7 = L L ^AL L d Rl ,d R 
{ ^ = LlYAL l u r1 ,u r 



Of L = (D^e R \ e AD,L L 
0%>=e R \ e AL L Q L \ D d R (4) 



(5)_7 , MAr ^ „ » n RL t'R'^e L ""R^U ' VL 



o£ = e^A e AL L « RCV A^T 2 g L 



and 



G>g = L L y AL L e^e* (5) 



0<5 = 5 nj CLirr a Li L™yV a L£ 



1 Throughout this paper we use four-component spinor fields, and ip c = —ij 2 tp* denotes the charge 
conjugate of the field ip. We also use v = (H°) ~ 174 GeV. 
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The MLFV hypothesis fixes the couplings A and 5 to be: 

f^r- Umf.U^ minimal field content 

A = { M " ~ ( 6 ) 

I Um v W extended field content, CP limit 

and 

rT I U*m u W minimal field content , , 

5 = 5 T =<Z~ ~+ (7) 
\^U*m v W extended field content 

in terms of the diagonal light neutrino mass matrix m v and the PMNS matrix U. To 
arrive at the above results we have used that \ e and assumed perturbative behavior 
in A and 5, stopping at leading order in (A LN m u )/v 2 (and (M u m u )/v 2 ). Therefore, for 
consistency, in the extended field case the operators should be dropped. 

The explicit form of the couplings Ay in terms of neutrino mass splitting and mixing 
angles (using the parameterization of the PMNS matrix U of Ref. |5J) is reported in 
Ref. PP for both minimal and extended field scenarios. Indicating by s and c the sine 
and cosine of the solar mixing angle and by S13 = sin #13 [S], the explicit form of Sij in 
the case of minimal field content (the only case we need) to first order in S13 is, 

At,N r 9 _„•„,. o 1 At 



<• * V LJN r 2 —ion 2 — iao 1 J1 LJN , 

See = —5- [c e 1 m Ul + s e 2 m U2 \ = d ee 



j- Aln 1 r -i a , -i ao is l Aln , 

Oeu, = —~ — 7= -see + see 2 m V2 + s 13 e m U3 \ = d eil 

v 2 V2 " v 

e Aln 1 r 2 -iai 2 -iao , 1 Aln , / Q \ 



Aln 1 rJ2 „_ iaim , „2„-ia 2 _ , _ 1 _ Aln 
v 2 2 

Aln 1 r - ia , -ia 2 is 1 A LN 

[see — see 2 m„ 2 + Si 3 e m U3 \ = d e 



v 2 y/2 1 

Aln 1 r 2 -iai 2 -iao 1 Aln , 

V = —5-7: I -s e 'm„, - c e 2 m^ 2 + = d Mr 



Throughout this paper we will assume maximal atmospheric mixing and the following 
central values for the remaining neutrino mixing parameters [S]: Am^ = 8.0 x 10~ 5 eV 2 , 
Am 2 atm = 2.5 x 10" 3 eV 2 , # sol = 33°. 



3 t — > I' I" I'" at tree level in the effective theory 

Within the effective theory described in the previous section we can now study decays of 
the type i — > £'£"£"'. At tree level these decays receive contributions from the transition 
magnetic moment operators (photon exchange), the four-lepton operators 0^1, 

and from (Z exchange). 



3.1 Rates 



The integrated decay rate for \i 
over final polarizations is: 



eee averaged over initial polarizations and summed 



r 



fi — >3e 



r 



t» 4 |A 



e/x I 



/i—*euu 



A 4 

LFV 



2|a_| 2 — 8Re(aQa_) — 4Re(a,Qa + ) + 6/|a | 



(9) 
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where Y ^ evV = m^/(t> 4 1536 7r 3 ), the a« coefficients are 
a + = sin 2 ^(ii + C S) + C g 

n l \lJX) i | Jl) | J2) 2( ^ee , J4) (5) \ /, n x 

a_ - (sin # w - -J(c iL + c LL ) + c AL + c 4L + — (c 4L + c 4L ) (1U) 



a = 2e 2 (4 1 | - cgi)* 
and / is a phase space integral given by (see appendix [X] for details) 

V^-|ln(2)-? + (^ . (11) 



3 \m e J 3 9 

Note that a_ depends not only on the Wilson coefficients but also on the ratio of effective 
FCNC couplings mediating fj, — > e transitions, namely S efl S2 e and A efl . We shall study 
the consequences of this in Section l3~3l 

The rate for r — > fifip, and for r — > eee is given similarly, replacing A eA1 — > A^, 

Vee V 5 w and Ae M V« -> S er S *ee, respectively, and m M -> m T . 

The integrated rate for t~ — > e~/i~/i + is given by 



r - = r 



v 4 |A er , 2 



a; 1 



a+| 2 + |a_| 2 - 4Re[a*(a+ + a_)] + 127|a | ( 12 ) 



LFV 

where a+, a are given in Eq. (JTUJl while 

5- = (sin 2 ^ - i)( C « + eg) + cffi + eg + ^ (eg + 45) , (13) 

and J is a phase space integral reported in appendix El The rate for r~ — > fi~e + e~ is 
obtained by exchanging the labels e <-> /i in A and 5. 

Finally, the rate for r~ — >■ ji~ n~e + is given by Eq. (jHJ) replacing \i — > r, but with 



a = 

( C 4L + C 4l!) 



^e.^, (4) , (5), (14) 



A e . 

(no photon and Z exchange contribution) and similarly for r _ — > e~e~ fi + , exchanging 
the labels e <-> /i in 5. 

In order to explore the phenomenology of these results, it is convenient to express 

the rates in terms of dimensionless couplings with the large scale Aln (for minimal 

field content) or M v (for extended field content) factored out. Hence we introduce the 
quantities a^, bij and dij defined through 

{A 2 
-jr-ciij minimal field content 
^bij extended field content 
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and Eq. (jHJ). These quantities depend only on low energy masses and mixing angles, 
and can be readily estimated. As an example we report the resulting expression for the 
[A — ► 3e rate: 

r M -*3e/r^ ei ,p = \a + \ 2 + 2|a_| 2 - 8Re(a* a_) - ARe(a* Q a + ) + 6/|a | 2 



x 



x 



'jkN_V| a 12 

Alfv/ 1 l" e ^l 



minimal field content (16) 



j^i) I Kfi\ 2 extended field content 



3.2 Phenomenology: extended field content 

Let us consider first the MLFV realization with extended field content, where the dis- 
cussion is somewhat simpler. In this case the contributions of Ojj^ and Ofl to fi — > eee 
are negligible, since they involve the combination S^ e 5* e which is suppressed by a small 
neutrino mass relative to the contributions to the amplitude that are proportional to 
A Me . As a consequence, the transition between families i and j is governed just by 6^-, 
the only allowed FCNC effective coupling. Moreover, for Wilson coefficients of 0(1), 
the expression in square brackets in Eq. (JTfij) is 0(1) and therefore, for a given value of 
{vM v j Alfv 2 ), the branching fraction for /j — > 3e is determined by \b e ^\ 2 . 

Fig. shows the behavior of \b eil \ 2 as a function of (the poorly determined mixing 
angle) S13 varying the lightest neutrino mass in the experimentally allowed range < 
^min 5; 0.2 eV jS], in the case of normal ordering of the neutrino mass spectrum. The 
plots correspond to the two CP conserving values of the phase 5 in the PMNS matrix: 
5 = (left panel) and 5 = n (right panel). As can be seen, the present uncertainties 
in s 13 and the absolute scale of neutrino spectrum induce a variation of \b e ^\ 2 over a 
couple of orders of magnitude. Similar results apply in the case of inverted spectrum. 
For M v = 6 x 10 7 A LFV /w, which saturates the perturbative Yukawa coupling bound 
on M v when Alfv = 50 TeV, one obtains \i — > 3e branching fractions of order 10~ 12 , 
comparable to the 90% CL. limit of 1.0 x 10~ 12 7J. The r — >• 3/i branching ratio is 
governed by IVI 2 , which does not depend on S13 and for given m min is typically two 
orders of magnitude larger than |6 e /J 2 , leading to B T ^ 3fl ~ 10~ 10 . This is below current 
experimental sensitivities jHl E] ■ 

The inclusion of 4-lepton processes in the phenomenological analysis offers more ways 
to test the hypothesis of minimal lepton flavor violation. We report below two noteworthy 
features, which with appropriate changes (see next subsection) extend to the minimal 
field scenario as well. 

• As pointed out in pQ , testable predictions of MLFV involve ratios of FCNC tran- 
sitions between two different families (e.g. /i — ► e vs r — > /x). In the case of 
three-lepton final states, several operators contribute and they pick a different 
phase space weight, in principle. So predictions are not as clean as in the case of 
lepton-photon final state. We find 





+3e 


\b ef x 


2 1 


^3e 


B r . 


^3/i 




2 1 


->3/x 



w ; (17) 



J T^3^ 
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Figure 1: \b e ^\ 2 as a function of S13 for normal neutrino mass hierarchy and the PMNS 
phase 5 = (left) or 5 = it (right). The shaded bands correspond to variations of the 
lightest neutrino mass in the range < m min < 0.2 eV. The upper edge of the bands 
corresponds to m m i n = 0. 



where the phase space integrals are i^-^e = 9.886, I T ^n — 3.26 and r c is the 
following combination of Wilson coefficients (in the absence of contact operators 
r c = 0): 

r c = —— (\a + \ 2 + 2\a„\ 2 - 8Re(a* a_) - 4Re(a*a+)) . (18) 

O Go 

Modulo the last (r c dependent) factor which is of 0(1), Eq. ([17]) has the same struc- 
ture of the prediction for the ratio of radiative decays B^ ei / B T ^^ n = |& eM | 2 /|6 MT | 2 , 
and will offer another way to test the MLFV pattern B T ^^ 3> B^ e ~ B T ^ e . A 
precise test requires knowledge of the ratio r c , which can be determined if additional 
LFV modes are observed, as we now discuss. 

Ratios of FCNC transitions between the same two families (e.g. \i — > ej vs \i — > 3e 
or r — > ^7 vs r — > 3/i and r — ► //ee) are determined by known phase space factors 
and ratios of various Wilson coefficients. If both fi — > and \i — > eee processes are 
observed one can begin to disentangle the effects of photon exchange from those 
of contact interactions, through the ratio of rates 



r^ 3e _ ct 

f ~ 4^ ^ 3e X 

1 fi—te-y 



1 + 



-'/li— >3e 



(19) 



with I^3 e defined in Eq. (jTTj) and r c in Eq. (JTHJ). The above ratio reduces to 
a/47r I^ze in the absence of contact interactions and offers a way to experimentally 
determine r c . One obtains the same result for r r _ >3M /r T _^ jU7 (with I (U _ i .3 e — > / T ^3/x 
and the same r c ) and a similar one for T T ^ iieE /T T ^ m . 



3.3 Phenomenology: minimal field content 

In the case of minimal field content, the analysis of 4-lepton processes is complicated by 
the contributions of and Of^ , implying that there are two FCNC effective couplings 
that mediate \i — > 3e: a Me and d^d*^ (see Eqs. (jOJ) and (jTHI )- This naturally leads us to 
consider two cases. 
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0.05 0.1 0.15 0.2 0.25 

Sl3 



Figure 2: Range of | function of S13. The shaded band corresponds to varying 

the PMNS phase S between and 2n. The range of \a eT \ 2 is the same as this, but is 
anti-correlated: the maximum of |a ep | 2 , achieved at 5 = 0, corresponds to the minimum 
of \a eT \ 2 . 



(i) If the dynamics of the underlying MLFV model results in c\£ 

4, 5, then the analysis of /1 — ► 3e processes parallels the one of the extended field content 
given in the previous section, with the changes \b^ e \ 2 — > \a^ e \ 2 and (wM^/Alfv 2 ) — > 
(A ln /A L fv) 2 - As before, for a given value of A L n/Alfv ; |<Ve| 2 determines the /i — > 3e 
branching ratio. We show in Fig. |21 the range of allowed values of | function of 

S13, varying the PMNS phase 5 between and 27r. We see that, except for very small 
S13, in the minimal field content scenario a ratio Aln/Alfv = 10 10 implies a branching 
fraction for fi —>■ eee as large as a few times 10~ 12 , comparable the 90% CL. limit 
of 1.0 x 10 _12 j5]. The range of \a eT \ 2 is the same as for |a eAt | 2 , but is anti-correlated: 
the maximum of |a eAt | 2 , achieved at 5 = 0, corresponds to the minimum of |a er | 2 . The 
quantity a Mr is independent of S13, la^rl 2 = 1-6 X 10 -51 . In this scenario Eq. (fTTj) remains 
valid with the substitution \hj\ 2 —> \ciij\ 2 , and r c (given in Eq. (fTKj) ) involves just ratios 
of Wilson coefficients that can be determined experimentally by considering the ratio 
r^3e/r M ^ e7 as in Eq. flH|. 

(ii) If the underlying dynamics is such that all Wilson coefficients are of comparable 
magnitude, then it is important to address the relative size of the FCNC couplings a^ e 
and d^ e d* ee , as there is a potential competition between them. In Fig. Elwe plot the ratio 
r efM = \2d e ^dlJ a efl \ as a function of m min for < m min < 0.2 eV, in the case of normal 
hierarchy (left) and inverted hierarchy (right). The points correspond to a scanning in 
parameter space with < S13 < 0.20 and < a\^, 5 < 2n for the CP violating phases in 
the PMNS matrix. For a given value of m m j n , the upper edge of the range corresponds 
to Si 3 — > 0. As can be seen, r^ e displays a very strong dependence on the absolute 
mass scale of the neutrino spectrum. Typically we find r efJL 3> 1, except in the normal 
hierarchy case and for m m i n — > 0. Moreover, for a given value of m min there are regions 
of parameter space where cancellations induced by CP violating phases drive r Me to zero. 

In the general case (ii), Eq. (fT9*j) remains valid but now r c of Eq. (fTSj) depends on 
ratios of Wilson coefficients and the ratio of FCNC effective couplings r eA1 = \2d e ^d* ee / a e ^\. 
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NORMAL HIERARCHY INVERTED HIERARCHY 




m m in (eV) m min (eV) 



Figure 3: Range of the ratio |2d eA1 <i* e /a eM | as a function of m min (eV) in the case of 
normal hierarchy (left) and inverted hierarchy (right). We use the range < S13 < 0.20 
and we allow the CP violating phases in the PMNS matrix to vary between and 2n. 
For a given value of m min , the upper edge of the range corresponds to S13 — > 0. 



The analogue of Eq. ()17j) now reads 



1 r c {r e/1 ) 

-B^—fSe I *^e/x I -^At— >3e -^/i— >3e 



(20) 



It- 



►3p 



If one could measure both r ^ 3e /r ^ n and r r ^3 M /T T ^ M7 (which determine r c (r eAi ) and 
r c( r fir) respectively), then Eq. (|2Uj) provides another way to test the MLFV hypothesis, 
insensitive to specific model details. 



4 Disentangling operators: hadronic processes 

As data becomes available on different lepton flavor violating processes it will be possible 
to begin disentangling the contributions of different operators. The rates for fi — > cy 
and —>■ 3e decays and for /x-to-e conversion in nuclei depend through different linear 
combinations on the Wilson coefficients in the effective lagrangian in Eq. Q. The 
ratio r M ^3 e /r M ^ e7 is sensitive to contact four-lepton operators already present in the 
operator basis (0^1) or induced by O^f via Z° exchange (see Eq. (fTH|) ). Similarly, 
r M ^ e conversion /r A1 ^ e7 is sensitive to the contact terms Of? L involving quarks either directly 
(Oll ^) or v ^ a Z° exchange (O^ ). In absence of contact operators the latter ratio is, 

^ = rrD A (21) 

where we use the notation of Ref. ^Hj for the dimensionless nucleus-dependent overlap 
integral D A . Typical values for the overlap integral in the case of a light and heavy 
nucleus are (TTF: D Ai ~ 0.036 and D Au ~ 0.19. A deviation from the result of Eq. (|2"Tf 
could be attributed to the contribution of Oyl. 
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C H -J 

0.05 0.1 0.15 0.2 

m min (eV) 

Figure 4: Ratio T(r — > nfi)/T(r — > jj,uu) in units of 10~ 27 (M u v / Al FY ) 2 , in the extended 
field content scenario as a function of the lightest neutrino mass (in eV). The red (lower) 
curve is for the normal hierarchy of neutrino masses while the green (upper) curve is for 
the inverted hierarchy. 



After the cancellation of the MECO ^T] experiment, which would have had a single 
event sensitivity to branching fraction for /i-to-e conversion of 2 x 10 -17 , it seems that 
experimental searches of fi-to-e conversion are in the far future. This motivates us to 
explore the sensitivity of other processes to Of^\ In this section we study other decays 
where hadrons are involved and contribute at tree level. Finally, in Sec. 03 we 

shall investigate whether can gi ye significant contributions to /i — > 3e through 

loop effects in models where the 4L operators are not generated at the scale Alfv (and 
therefore their coefficients can be discarded). 



4.1 r 



7T°/i 



Consider r — > fiX, where X = ir , p , etc. With a non-vanishing left handed neutrino 
Majorana mass term, electroweak interactions will produce a non-zero amplitude at one 
loop. Roughly, this amplitude corresponds to an operator 

1 {Uml&)^ L YT L ){q la { lh )q). (22) 



O 



SM 



a 



s 2 M 2 M 2 

b w lvl W lvl Z J 



This is much smaller than the contribution from the operators generated at the LFV 
scale. For example, for minimal content the factor in square brackets, of order 10~ 10 GeV -4 , 
is replaced by 1/f 4 (Aln/Alfv) 2 ~ 10~ 9 GeV~ 4 (ALN/ALFv) 2 , which for any reasonable 
value of the ratio Aln/Alfv is orders of magnitude bigger. 
A straightforward calculation gives 

' m l^lrfn r (1) (2) _ (Ad) (Au) „ (5) 1 2 
' [ C LL i* C LL C LL i* C LL ' AC LL\ > 



/17T 



2567T 



A 4 

LFV 



(23) 



where f n = 130 MeV is the ir decay constant. Numerically we find, in the minimal model 



fxir 



T(r — > \ivv) 



5.2 x 10~ 52 



A 



LN 



A 



LFV 



LL 



+ C 



(2) 
LL 



C 



(Ad) (Au) 



LL 



+ C 



] 2 



(24) 
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Using Aln/Alfv = 10 9 which corresponds to Br(/i — > ej) ~ 10 13 , and assuming the 



coefficients c^ L are of order unity, this yields Br(r — > jm ) ~ 10" 



15 



For the extended model we obtain 



T(r -> /ivr°) 
r(r — > |Uz/i/) 



5.0 x ID" 27 ' vM " 



A 2 

LFV 



(1) (2) _ (4d) (4u) 
LL "r LL LL LL 



(5)1 2 



(25) 



v/A 



m* tm . The result is 



in the normal hierarchy case, assuming m min = m Ul <C m U3 
plotted as a function of the lightest neutrino mass, m m i n , in Fig. 0]both for normal and 
inverted hierarchies. If M u = 3 x 10 5 A 2 FV /t>, which again corresponds to Br(/i — > ej) 
1Q- 13 , we obtain Br(r ^tt°) ~ 10~ 15 . 



4.2 7T° -> /i+e" 

With almost no extra work we can compute 



r(7r° -> /x+e" 



f* m yi( m *- m ») l A e M | r (1) (2) _ (4d) (4u) , Q (5) 1 2 

128vrm3 A4 py L c ll + c iL c iL + c LL + zc LL \ 



Numerically, using l/r(ir°) = 7.83 eV, 



(26) 



Br(.° p+e-) = 1.4xlO-»[c<'» + c« 2 »-cif + c«t'+2cS] 2 1 TTh 

[(M^/A LFV ) 2 |6 eM | 2 extended 

(27) 



Assuming ~ 0(1) and Aln/Alfv = 10 9 one expects Br(7r° —> fi + e 
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4.3 J lij) -> r/i and T 



The lepton flavor violating decays of heavier neutral mesons are interesting because 
decays into r's become energetically allowed. By normalizing the decay rate for V — > ££' 
(for V = J/i/j, T) to the well measured — > ee we can make a prediction free of hadronic 
uncertainties: 



r(v 



ee 



e" 1 V 2v / 



(Aln/Alfv) 4 !^. 
(M^/A LFV ) 2 |6, 



12 
12 



with 



v 




+ c (2) Hi - - sin 2 1 + c (3) + c (4n) - c 

^ C LlH2 3 b111 C LL ^ C LL C 



(3) 



„(4tt) 



.(5) 



LL 



+ C 



(2) 
LL 



(-i + |sin 2 



W ) + 



S3) 
'LL 



+ C LL + C 



(5) 
LL 



minimal 

extended 

V = J/tj) 
V = T 



In particular, for T decays into r\x this gives 



Br(T 



.(1) 
'LL 



+ C 



(2) 



LL 



1 2. 2/1 
-2 + 3 am 



.(3) 



(4d) 



+ C LL + C >.r. + C 



LL 



,(5) 
-LL 



X 



X 



3.2 x 10^ 55 (Aln/Alfv) minimal 

3.2 x 10" 30 (M„v/A LFV ) 2 extended, normal h, m un 

4.3 x 10~ 30 (M l/ t)/A 2 jFV ) 2 extended, inverted h, m h 



(28) 



(29) 



(30) 



0, (31) 
= 0. 



11 



Assuming again c[ L ~ 0(1) and Aln/Alfv = 10 9 one expects Br(T — > t + n~) ~ 1CT 20 . 

All the hadronic modes we have considered here, for acceptable values of the ratios 
Aln/Alfv and vM^/Alfv 2 and of the Wilson coefficients, are predicted to be well below 
foreseeable experimental sensitivities |12j . On one hand this result can be considered 
as a generic prediction of minimal lepton flavor violation and can be used in the future 
to test this framework (e.g. observation of T — ► t + jjT at B factories would strongly 
disfavor it). On the other hand we see that within minimal lepton flavor violation the 
most sensitive probes of are: (i) fi — > e conversion in nuclei and (ii) purely leptonic 

processes such as // — > 3e (to which O^ 5 ^ contribute via loops). 

5 \i —> eee in the absence of 4L operators 

The operators can give rise to /j, — > eee via one loop graphs. These contributions 

can be significant in models where the 4L operators are not generated at the scale Alfv> 
or are generated with tiny Wilson coefficients. We will analyze now this dynamical 
scenario, which amounts to the assumption c^(A LFV ) <C c21(A L fv)- in order to address 
the sensitivity of /i — > eee to , we organize our analysis as follows: first we use the 
renormalization group to evolve the effective lagrangian from /j ~ Alfv down to fi ~ 1 
GeV. Then we take /i — > 3e matrix elements of the relevant operators, using perturbative 
QCD to deal with the heavy quarks and chiral perturbation theory to deal with the light 
quark loops. 

5.1 RG analysis 

We use the renormalization group equations (RGEs) to determine the low energy effective 
lagrangian. This is done is three steps: first, the RGE in the unbroken phase of the 
SU(2) x U(l) theory is used to compute the coefficients in the effective lagrangian in 
Eq. (J3J) down to a scale \i ~ M z . In the second step the coefficients are matched to those 
of an effective lagrangian for the theory in the broken symmetry phase of SU(2) x U(l). 
And third, the coefficients of this effective lagrangian are computed at \i ~ 1 GeV using 
the RGE for the theory with only U(l) gauge group. The calculation is rather lengthy 
and technical, so the details are presented in Appendix IB1 

In the first step of the calculation we compute for M z < fi < A L fv the mixing of 
coefficients c^ L into the coefficients of the four-lepton contact interactions, cfl, assuming 
Ql(Alfv) *C cll(Alfv)- In the second step, an effective theory with U(l) gauge field is 
constructed with effective lagrangian given by matching at /i ~ Mz to the theory with 
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SU(2) x U(l) symmetry. The relevant operators in the new effective theory are 





o 



o 



(3) 
LL 

{id) 
LL 

(4») 
LL 



o 
6 



(5) 
LL 

(1) 
RL 



o 



(3) 
1L 



e L ^Ae L (ulI^ul + d L ^^d L ) 

e L j^Ae L d R j^d R 

e L ^Ae L ur^^ur 

-e L ^Ae L (ulj^ul - d L j^d L ) 

eve R a^X e Ae L 

e L j^Ae L tLl^L 



(32) 



where sums over quark generations are understood (the top quark has been integrated 
out), and the effective lagrangian is 



C 



A 2 

LFV 



s(0 Ad) I ^W/SW 

C LL W LL 4 T ■ 



+ 



A 2 
LFV 



E 



416g + h. c . 



(33) 



The coefficients cll,rl,al , obtained by matching, are given in Eq. (jo3j) of Appendix |Bj 
Finally, in the third step we solve the RGE satisfied by c^l and cfl- Ignoring the 
running of cll/rl (that is working to order g 2 ), we finally obtain the following low scale 
(/i < M z ) Wilson coefficients: 



4V) = - 



S (3) 



In 



,(3) 
-LL 



+ 4cS ) -3cVT-7c 



(4d) 



+ 



.(3) 



,(4u) 



LL 



.(5)' 
-LL 



+ Cr r/ - C 



(A*) = -r 




(5)' 
LL 



+ 2c 



(4m) _ (4d)\ 



3c 



LL 



(4d) 
LL 



-LL 



7r: 



(5) 
LL 



+ C> r/ - C 



(5) 



LL 



+ 2c 



(4u) _ (4d)\ 
LL C LL J 



/ 2 sin 2 9,, 



.(5) 
-LL 



|+sin 2 ^ (eg] 



(!) , J 1 ) 



+ C 



LL 



(34) 



5.2 Matrix elements 



At the order in g l we are working, the matrix elements of 0\£ have to be taken at tree 
level, while those of L \ have to be evaluated at one loop. Consider for definiteness the 
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e e 



Figure 5: One-loop contribution to (3e|Oj^J//). The heavy quark contribution is treated 
in perturbation theory, while the light quark contribution is evaluated using chiral per- 
turbation theory. 



decay fi — > e (pi) e (p 2 ) e + {p E ). Let us define: 

M^(pi,p 2 ) = A £M m(pi)7 m P l m(p m ) • u(p 2 )iuv{pe) M(pi,p 2 ) = g^M^ipiipi) 

and the kinematic variables x = mf 3 /m 2 = (p ei + Pe) 2 / 772 ^ and y = m^/m 2 ^ 
(Pe 2 +Ps) 2 /m 2 l . 

Matrix elements of 0^ L : heavy quark contributions 

Considering the graphs of Fig. |2]with heavy quark internal loops, we find 



<3e|VcM3>> 



a 



71 



(2F c {y)-F b {y))cf L + 2F c {y)cT L 



,(4«) 



- F 6 (y)6gf - (2F c (y) + F b (y))cf L ] x M( Pl ,p 2 ) ( 35 ) 



Pi <-> P2 ; y 



X 



} 



where F q (x) = F(m q , {m q /m^) 2 fx), with 



F(m, z) 



' | In 2 -f - ± _ | + i (1 + 2.) arctan 



In 



1—^/1=42 



+ i7T 



In muon decays m 2 3 , mf 2 ~ g 2 < m 2 , and by using 



for Az > 1 

for 4z < 1 
(36) 



F(m, m 2 /q 2 



1 , 2m 2 

-log— 2" 



30m< 



+ 



the heavy quark contribution can be safely approximated to a constant (local) term. 
Using the q 2 /m 2 — > limit and the notation F g = 1/3 log [(\/2ra g )///] we then obtain: 



7T 



X 



(2F C - F b )cf L + 2FJ%$ - F b c^ - (2F C + F b )cf L 
+ (Of2) 



(37) 
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It is easy to check that the scale dependence of the above matrix elements cancels the 
one induced by b, c loops in the Wilson coefficients c^, cfl given in Eq. 

Matrix elements of 0^ L : light quark contributions 

The operators 0^ L have the structure 

1 

2 

where the flavor structure of the currents Vp , can be read off Eq. Only the 

vector current contributes to the matrix elements in question and we find: 

.2 1 



6f L = e L ^Ae L £ 



yii) ± ^(0 



(38) 



^e\of L \^) = -C-M^{p u p 2 ) 



(p 2 + Pe 



■X 



x / d 4 x^ P2+p ^(0\T[V;{0)V^ M (x))\0) - {pi^p 2 } 



Defining then 



T (V*(0)Vr(x) |0> 



(4tt)< 



( M 2 " q*v) ■ ^' EM (Q 2 ' 



one arrives at: 



(39) 



(40) 



<3e|0S» 



2(4vr) ; 



(41) 



In order to evaluate the VV correlator at low momentum transfer we use £77(3) x SU (3) 
chiral perturbation theory to order p 4 . Evaluating the one-loop diagrams depicted 
in Fig. |H1 (with vertices from the 0(p 2 ) chiral lagrangian) and adding the counterterm 
contributions from the local 0(p A ) effective lagrangian ^Hj, we find: 



n 



4u,EM 



i=iv,K 



f(z,m 2 /f£ 



_I l n Il£ + ^ _ » + 2 (4 _ 1)3/2 cot -i(V4^T) . 
3 u£ 9 3 3 V ; V ; 



(42) 



The effective couplings L\o and H\ are defined in Ref. and they cancel the loop- 
induced dependence on the chiral renormalization scale \x x . The constant H\, however, 
depends on the renormalization scheme adopted in the short distance theory (at the 
quark level), which is in this case the MS scheme. This dependence is such as to cancel 
the \i dependence of the amplitude induced by the Wilson coefficients (see appendixlHlfor 
details). Similarly we find IT 3 ' EM (q 2 ) = 0, IT 4d ' EM (q 2 ) = -U iu ' EM (q 2 ), and IT 5 ' EM (g 2 ) = 
2II 4u ' EM (q 2 ). The final result for the light-quark loop matrix elements is 



a 1 
4^ 2 



C LL C LL 



w> + 2a ( 5 ) 



X 



x M( Pl ,p 2 )U^ EM ((p 2 +p s ) 2 ) { Pl ^ P2 } 



(43) 
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■ K 



7T, K 



V' 



EM 



(xywv 



Figure 6: One-loop contributions to the VV correlator denned in Eq. (|4"0J) in SU (3) x 
SU(3) chiral perturbation theory. The vertices come from the lowest order (0(p 2 )) chiral 
lagrangian coupled to external sources. 



Over the physical region for \i — > 3e decay the function n 4u ' EM (q 2 ) varies by less than 
1% 2 . Approximating it to a constant we find: 



(4-u) ~(4d) 



LL 



LL 



+ 2c 



(5) 
LL 



n 



4m,EM 



K) x (6«) + (6S> . (44) 



(3)\ 



In summary, the matrix element calculation results in: 

(3e| 6g>> = ( Khcavy _ q + K light _ q ) x ((62) + (62) 



(45) 



with 



^heavy— q 
^light— q 



a 

71 



m - F h )cf L + 2F c cT L > ~ F b ct a L > ~ (2F C + F h )cf L 



(4d) 



i(5) 



"8^ 



(4m) -(4<Z) 



LL 



LL 



+ 2c 



(5) 
LL 



n 



4u,EM 



K) 



(46) 



5.3 Rate 

The result in Eq. ()45|) implies that in the scenario considered in this section the /J — * 
rate is given in closed form by Eq. ©, replacing the coefficients a± of Eq. (flTjj) by: 



eee 



~^4L ^heavy— q 



-r {1) + 

u 4£ rb heavy- 



+ K light-q 
q "I - ^light— q 



(47) 



with given by Eq. (J54*j) . Taking into account the /i dependence of if[(/i x ,/i) it is 
straightforward to verify that a± do not depend on the renormalization scale [i. 

While the formulas given above are quite general, in order to illustrate at which level 
[i — > eee decays probe Of^\ let us assume that only one operator dominates. Moreover, 
let us focus on the operator contributing with the largest numerical coefficient, which 
turns out to be O ll . We find then: 



6-10- a x |<4 5 !(A LFV ) 




minimal field content 



b ef j\ 2 extended field content 



(48) 

2 For the low energy constants we use L r w {^i x = m p ) — —5.5 x 10 -3 (experiment, see e.g. |T3]) and 
|iJ[(/i x = m p , /J, — lGeV)| < l/(647r 2 ) (naive dimensional analysis JH| bound). 
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A comparison of Eqs. and (|16|) shows that the rates of I — > £'£"£"' decays in the 

scenario c^(Alfv) ^ c l1(^-lfv) are typically suppressed relative to their tree level 
counterparts by ~ 1CT 3 . This factor is larger than the naive estimate (a/ir) 2 because of 
the presence of sizeable logarithms. Similar logarithmic enhancements were first pointed 
out within a non-MLFV effective theory approach to lepton flavor violation in Ref. jTH] 
(in that reference the main focus was on the one-loop contributions to /i-to-e conversion 
in nuclei) . Finally, let us notice that apart from the overall suppression, in this dynamical 
scenario we still find the pattern T(t — > //) ^> T(fi — > e) ~ T(r — > e), which is dictated 
by the MLFV hypothesis. 

6 Conclusions 

The MLFV hypothesis provides a framework for discussing and analyzing the predic- 
tions of a large class of models where Lepton Flavor Violation arises solely from lepton 
mass matrices. The framework is attractive because it is both very general and fairly 
predictive. We have extended the results of Ref. [I], which introduced the MLFV hy- 
pothesis and analyzed /i-to-e conversion in Nuclei and radiative decays, e.g., \i — > cy, to 
decays involving hadrons, e.g., r — > fin, and to decays involving only charged leptons, 
e.g., ji — > 3e. To this end we extended the operator basis for the effective lagrangian 
introduced in pQ. In particular, we found five new purely leptonic operators, given in 
Eq. (JHJ), which were omitted in pQ because they do not contribute directly to the processes 
considered there. 

In the event that several LFV processes are observed and their rates measured one 
could begin to test the MLFV hypothesis since it predicts definite patterns of relative 
magnitudes of rates. We have therefore computed the rates for the hadronic processes 
t — > Etc {I = fi, e), tt° — ► /i + e - and V — > T/i (V = J/ip,T), and for purely charged 
lepton decays // -> 3e, r -> U (£ = e,/x), r -> £££' and r -»■ £££' ((£,£') = (e,/x), (/i,e)). 
One definite prediction of the MLFV hypothesis is that the rates for decays involving 
hadrons are exceedingly small. For example, if the scales of LFV and LNV are such 
that Br(/i -> ey) ~ 10" 13 , then Br(r fin ) ~ 10" 15 , Br(7r° -> n+e~) ~ 10^ 25 and 
Br(T — > r + fi~) ~ 10 -20 , well below the sensitivity of foreseeable experiments. 

We analyzed the decays involving only charged leptons first at tree level and then 
through loop processes, in order to investigate the contribution of the operators involving 
quarks (under the assumptions that the coefficients of the 4-lepton operators of Eq. ((SJ) 
are negligibly small). Some general conclusions apply to both cases: (i) LFV 3£ decays 
display the same pattern B T ^^ ^> B^_ e ~ B T ^ e as the radiative decays, and hence offer 
an alternative way to test the MLFV hypothesis; (ii) ratios of transitions between the 
same two families, as, for example r(/i — > 3e)/r(/i — > ej), are determined by ratios of 
various Wilson coefficients and therefore the combined measurement would be valuable 
in disentangling the contributions of different operators. 

More specifically, the tree level rates for decays involving only charged leptons can 
be very significant. If the coefficients in the effective lagrangian are of order 1 and the 
ratio of scales of LNV and LFV is as large as allowed by the condition that the Yukawa 
couplings are perturbative, then in the case of minimal field content Br(fi — > 3e) could 
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easily exceed the present 90% C.L. limit of 1.0 x 10 -12 . In the extended field content 
case, under the same assumptions, this branching fraction is comparable to the current 
experimental limit, and the fractions for LFV r decays are three orders of magnitude 
smaller than the current limits of 1 — 3 x 10 [HI 10]- Note that while the rates for LFV 
t decays are larger than for fi — > 3e, the higher experimental sensitivity of the latter is 
on the verge of placing stringent bounds on the theory and therefore it would be highly 
desirable to pursue higher sensitivity in this mode. 

If the new physics giving rise to the effective lagrangian at the LFV scale is tied to 
the quark sector in such a way that only operators involving quarks and leptons together 
were produced by this dynamics, then the decays to 3£ would not proceed at tree level. 
To investigate this scenario we computed the rates for 3£ decays through loops, assuming 
the coefficients of the 4-lepton operators of Eq. (jSJ) are negligibly small. Alternatively, 
since the rates for processes involving hadrons are exceedingly small, a more sensitive 
probe of them may arise through their loop effect on 3£ decays. Not surprisingly we 
found the rates are suppressed relative to their tree level counterparts by ~ 10 -3 . The 
calculation is interesting in its own right. It involves the computation of a low energy 
effective lagrangian obtained by integrating out the heavy quarks and the W ± and Z° 
vector bosons, and also requires the use of the chiral lagrangian to properly describe 
the low energy physics involving light quarks. As a side result a new sum rule for the 
Gasser-Leutwyler counter-term H{ was derived in Appendix O 

A discrimination of various dynamical scenarios within MLFV will be possible only 
when more than one LFV decay mode is observed. This remains true even when analyz- 
ing data beyond the minimal flavor violation hypothesis [TZj. We therefore emphasize 
that it is highly desirable to complement existing experimental efforts jTHJ |Hl E] and 
pursue experimental searches of all the LFV /i and r decay modes. 
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A Differential distributions and phase space for £ — > 

eti" 

The differential rate for fi — > eee in the variables x = m^/m 2 ^ = (p ei + pi) 2 jm 2 ^ an d 
y = m 2 3 /m 2 = (p e2 +Pe) 2 /m 2 is: 



d 2 T m 5 1 

[2 (2x 3 - 2x 2 + x + 2y 3 - 2y 2 + y) \a \ 



dxdy 2 9 7t 3 Al FV xy 

+ |a_| 2 (-8yx 3 - 16y 2 x 2 + 8yx 2 - 8y 3 x + 8y 2 x) + \a + \ 2 (-2yx 3 + 2yx 2 - 2y 3 x + 2y 2 x) 
+ a ((8yx 2 + 8y 2 x - 8yx) a*_ + (-2yx 2 - 2y 2 x) a* + ) 

+ (a+ (-2yx 2 - 2y 2 x) + a_ (8yx 2 + 8y 2 x - 8yx)) a* ] (49) 

The phase space integral / appearing in Eq. |U] reads: 

/= / dx dy{ — [y(l-x)+x(l-y)-2(l-x-y)(x 2 + y 2 )]+2(l-x-y)\ 
Jym in l x y J 

(50) 

in terms of x = (p ei + pi) 2 /m 2 ^ y = (p e2 +pe) 2 /m 2 l , and m = rue/m^. The integration 
limits are given by: 



1 

2/min — ~ 



1 — x + 3m =F 



2 __ i \X — 4m 2 ) [(1 — x — m 2 ) 2 — Am 2 x\ 



x 



(51) 



Numerically, J = 9.886 for /i — > eee I = 17.4 for r — ► eee and / = 3.26 for r — > fififi. An 
exact analytic expression for / cannot be readily obtained, but by expanding in powers 
of m the integrand in (j5Up after performing the y-integral, keeping the full dependence 
on m in the limits of x-integration gives the result in Eq. (fTTj) . which numerically gives 
10.4 for fi — > eee. 

The differential rate for r _ — > e~ fi~fi + in the variables x = m 2 2 /m 2 = {p e +Pfi-) 2 /ntf- 
and y = m 2 3 /m 2 = (p^- +p M +) 2 /m 2 is 



d 2 T ml 1 



[2 (2a; 2 + 2yx - 2x - y + l) \a \ 2 + (2xy - 2x 2 y) 



|2 



dxdy 2 8 7r 3 A£ FV y 

+ (-2y 3 - Axy 2 + 2y 2 - 2x 2 y + 2xy) \a + \ 2 + a ((2y 2 + 2xy - 2y) a\ - 2xya*_) 

+ (o+ (2y 2 + 2xy - 2y) - 2a^xy) a* ] (52) 

The integral over phase space /, appearing in Eq. (|12p. is given by (m = m^/m T ) 

. = ^ ^ l- y -2,(l-x- y) (53) 



The numerical value is 1.50 for r~ — > e~fi~fi + and 8.49 for the corresponding integral in 
the case of r~ — > fi~e~e + . Note that the limits of integration are different in these cases 
than in /x — > eee. 
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B Low Energy Effective Lagrangian 



In this appendix we give the details of the calculation described in Sec. EH that is, the use 
of the renormalization group equations to determine the low energy effective lagrangian. 
As explained there, this is done is three steps: first, the RGE in the unbroken phase of 
the SU (2) x U (1) theory is used to compute the coefficients in the effective lagrangian in 
Eq. (jSJ) down to a scale // ~ M Z - In the second step the coefficients are matched to those 
of an effective lagrangian for the theory in the broken symmetry phase of SU(2) x U(l). 
And third, the coefficients of this effective lagrangian are computed at \i ~ 1 GeV using 
the RGE for the theory with only U(l) gauge group. 

In the region fi ^ Mz (neglecting SU(2) x U(l) symmetry breaking) a one loop 
computation gives 

rfC 2 _ 1 AT AT Si r j3) , „(4«) J4d), 

LL \ 



,,ZZM±. - -N N — [r l ; 4- 2r ( ' - r ( 

dfi ~ 18^ cI ^4^ [Cll + ZCll Cl 

dc^ _ 1 ^_ (5) ( 54 ) 

a ^SL - -In N ^Mc (3) + 2c (4n) - c (4d) l 

where we have denoted by iV c the number of colors and by N g the number of quark 
generations. On the right hand side of these equations we have neglected terms propor- 
tional to c^l because we are considering a situation with C4l(Alfv) <S cll(Alfv)- A full 
solution of the RGE requires equations for c^ L (ji) too; however, since are small, we 
retain only up to order gf 2 in the amplitudes, and these terms drop out. 

In the second step, an effective theory with U(l) gauge symmetry is constructed with 
effective lagrangian given by matching at fi ~ Mz to the theory with SU(2) x U(l). In 
this process W ± , Z°, and top quark are integrated out. The relevant operators are given 
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in Eq. (|32|) . and the effective lagrangian in Eq. (}3*3*j) . Matching coefficients gives 
cf L {M z ) = cf L {M z ) - \ sin 2 e w [ C V(M z ) + cf L {M z )} 

+ \^QtF t N c [cf L {M z ) + c^{M z ) - cg(Mz)] 
^(Mz) = c£>(M*) - | sin 2 ^[cgCMz) + cg(M z )] 
+ |-g t F,iV c [ c ( 3 ) (M Z ) + c^\M z ) - cf L {M z )} 

A 71 

c ( ^(M z ) = cf L \M z ) + \ sin 2 6Ucg(M z ) + cg» (M z )] 

- ~g t F t N e [cg(Mz) + cg^M*) - cg(Mz)] 

3 7T / , 

1 1 (55) 
cf L {M z ) = cg(Mz) - (- - - sin 2 6y [c« (M z ) + cg> (M z )] 

- ~g*F t JV e [cg(M z ) + cfl\M z ) - cf L {M z )} 

I 71 

^l(Mz) = cg(Mz) + cg(Mz) - (i - sin 2 ^)[c«(Mz) + cf L (M z )} 

- -g t F,iV c [ci 3 i(Mz) + 4 u) (Mz) - cg(Mz)] 

7T 

\_„(3)/ il , x . .2 r (1)/,, x . (2), 



fiJ2(Mz) = cJS(Mz) + sin 2 ^[^(Mz) + c&(M z )] 

- ^F,iVjcg> (M z ) + cif(Mz) - cg(Mz)] 

4 1 i(^) = 4 1 i(^)-4 2 i(^) 

Finally, for the third step we need the RGE satisfied by cfl and cfl'. 



^ = ^dcfl _ e 2 iV ( 



dfi ^ d^i 9 ■ 4tt 2 



(n \-(3) , o "( 4u ) "( 4d ) i \-(5) 

(2n u - n d )4£ + 2n u c y LL ' - n d c LL ' - {2n u + n d )& L ' L 

(56) 

Here n u and n d again stand for the number of up-type and down-type quarks, re- 
spectively. It is easy to check that if one adds to the one loop amplitude from L l L , 
i — 3, . . . ,5, the tree level amplitude from 0\l , and one does not re-sum leading logs 
{i.e., one works to order e 2 only), then the \x dependence cancels in the amplitude (the 
/i dependence in c\£ cancels that of the loop amplitude). 
The solution to these equation is given in Sec. 15.11 in Eq. 



C A sum rule for H[(fj, x , fi) 

Let us return to the issue of matching the quark-level description with the chiral effective 
theory calculation. In particular, in this appendix we will show the explicit dependence 
of the chiral coupling Hi on the MS renormalization scale fi. 

The light quark contribution to the matrix element (3e|0^||/i) in eq. (jU]) is UV 
divergent, being proportional to U l,EM (q 2 ). In dimensional regularization the divergent 



21 



correlator reads (work with % — 4u for definiteness) 



n 



4u,EM/2 



NcQu 

(4tt) 2 



4 



L3 



- 7 + log 47T - 8F(m u , m 2 u /q 2 



(57) 



with F(m u , rn^/q 2 ) defined in Eq. (pjj) . The MS prescription amounts to subtracting the 
term proportional to — 7 + log47r). Writing II(g 2 ) = n^) + Tl(q 2 ) — n^), with 
q 2 in the regime of validity of perturbation theory, the MS subtracted correlator can be 
written as: 



(58) 
(59) 



-|-r4u,EM/ 2\ 4 N C Q U 
MS W ) 3 ( 47r ) 2 


log 


-H 


5" 
~ 3_ 


+ n 4 "' EM (g 2 ) - n 4M ' EM (g 2 ) 


4iV c g u 


log 


-2% 2 
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where we have used the notation fl(g 2 ) = n(g 2 ) — n(0). The first term has the correct 
short distance \x dependence. The remaining terms can be expressed through a (once 
subtracted) dispersion relation as 
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(60) 



If q 2 AySBi we can saturate the spectral function with the tttt and KK threshold 
contributions, up to higher order terms in 5 ,2 /A 2 5B . Doing so and comparing with the 
chiral perturbation theory result in eq. (|42j) leads to a matching condition involving the 
low energy constant L w and the contact term H\\ 
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A few comments are in order: 

• The LHS does not depend on fi x (due to chiral RG for Hi and L w ). 



+n 



4u,EM/„2 



(61) 



The RHS does not depend on q 2 . This can be seen by separating the dispersion 
integral for Il 4n ' EM (g 2 ) in an IR and UV region, and then using the free quark 
spectral function in the UV regime. This gives back a logarithmic dependence on 
q 2 that cancels the one in the other term. 

Eq. (joTj) explicitly shows how Hi depends on the short distance renormaliza- 
tion conventions and in particular on the scale /1. It also displays a purely non- 
perturbative contribution to the matching, namely n 4n,EM (gQ). In absence of reli- 
able esimates of n 4u,EM (g 2 ), in our numerics we use the experimental value for L 



r 
10 



and a bound for Hi{ji y 



m 



pi 



jj, = 1 GeV) based on naive dimensional analysis. 
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